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Abstract. Data obtained with the NICMOS instrument on board the Hubble Space Telescope (HST) have been 
used to determine the H-band luminosity function (LF) and mass function (MF) of three stellar fields in the 
^ ' globular cluster M 15, located ~ 7' from the cluster centre. The data confirm that the cluster MF has a character- 

istic mass of ~ 0.3 M©, as obtained by Paresce & De Marchi (2000) for a stellar field at 4f6 from the centre. By 
lO ■ combining the present data with those published by other authors for various radial distances (near the centre, 

at 20" and at 4'6), we have studied the radial variation of the LF due to the effects of mass segregation and 
derived the global mass function (GMF) using the Michie-King approach. The model that simultaneously best 
fits the LF at various locations, the surface brightness profile and the velocity dispersion profile suggests that the 
GMF should resemble a segmented power-law with the following indices: x ~ 0.8 for stars more massive than 
0.8 M , x ~ 0.9 for 0.3 - 0.8 M and x ~ -2.2 at smaller masses (Salpeter's IMF would have x = 1.35). The best 
' fitting model also suggests that the cluster mass is ~ 5.4 x 10 5 Mq and that the mass-to-light ratio is on average 

Oh' M/L v ~ 2.1, with M/L v ~ 3.7 in the core. A large amount of mass (~ 44%) is found in the cluster core in the 

form of stellar heavy remnants, which may be sufficient to explain the mass segregation in M15 without invoking 

'_h ' the presence of an intermediate-mass black hole. 
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1. Introduction tion. The observed stellar density profile can in fact be 



3 i reproduced by assuming the existence of either a central 

M 15 (NGC 7078) is a classical target for studying the in- h]a± hole (Guhathakurta et al . 1996) or a compact core 

ternal dynamics of globular clusters from the observed ag thfi byproduct of the cluster core collapse in the pres _ 

surface brightness and velocity dispersion profiles for it ence Qf diffuge dark m&tter (L&uer ^ &1 1QQ1) Vgry rg _ 

is at a relatively large distance from the Galactic plane Baumgardt et aL (2002) and Gerssen et aL (20 02) 

(Z G = -4.7 kpc; Harris 1996). This location, combined haye interpreted the kinematica i data obtained with STIS 

with an orbit of small ellipticity (e = 0.34; Dinescu et al. , w^nnr, j . . . • f . i 

* . J v '. and WFPC2 as due to either strong segregation of stel- 

1999), minimises the cluster interaction with the Galaxy , ± ij 

n J lar remnants (white dwarts and neutron stars, as already 

and hence tidal stripping and evaporation of stars from suggested by niingworth & King 1977) or to the presence 

the cluster outskirts. Ground-based observations revealed of & „ 1q3 m black hde in the core of M lg< Ag bef 



the presence of a central cusp which has been attributed both explanations are statist ically equivalent. 

to core collapse (Djorgovski & King, 1986). Recent HST j-, r ,. , , 

. i , , .„ , r further evidence of mass segregation m the central re- 

observations have not been able to clarify the nature of r , , ir r , , ,. r n ,. 

, . . , gions of M15 comes from the detection of colour gradi- 

the M15 cusp, m spite of their enhanced spatial resolu- , , , / TT D > , fT) , n , , , , , , 

r ' r ents whereby (U-B) and (B-V) colours get bluer towards 
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Table 1. Log of the observations 





Dataset 


RA (h) 


DEC (d) 


PA (d) of the 


Filter 


Number of 


Total exposure 






J=2000 


J=2000 


detector Y axis 




images 


time (s) 


Field 1 


n4k6r0hca 


21:30:17.22 


12:15:55.4 


0.341 


F110W 


4 


767.835 




n4k6r0hda 


21:30:17.22 


12:15:55.4 


0.341 


F160W 


4 


1023.823 


Field 2 


n4k6r2hsa 


21:30:19.01 


12:15:51.5 


0.342 


F110W 


14 


2687.426 




n4k6r2hta 


21:30:19.01 


12:15:51.5 


0.342 


F160W 


14 


3583.384 


Field 3 


n4k6u0v8a 


21:30:15.10 


12:16:34.2 


-8.546 


F110W 


18 


3455.262 




n4k6u0v9a 


21:30:15.10 


12:16:34.2 


-8.546 


F160W 


18 


4607.208 



driving mechanism. Dc Marchi & Paresce (1994) have re- 
solved with HST/FOC a large number of bright blue stars 
in the core of M 15, the majority of which can be classi- 
fied as blue stragglers. Nevertheless, the rest appear to 
belong to a new, as yet unidentified class of very blue 
stars. Amongst several possibilities (such as Early-Post 
AGB, subdwarfs and well- mixed single stars), De Marchi 
& Paresce (1994; 1996) have suggested that dynamical in- 
teractions and close encounters could have stripped off the 
envelope of red giant stars, enhancing their mass loss and 
heading their evolution towards the late stage of helium 
white dwarf and CO white dwarf. 

Because of minimal interactions with the Galaxy, the 
outskirts of M 15 have been observed to constrain the clus- 
ter initial MF (IMF) of low mass stars. For example, De 
Marchi & Paresce (1995) performed deep HST/WFPC2 
photometry of a field 4.'6 NW of the centre and derived 
the LF of main-sequence stars down to Mj ~ 10, or two 
magnitudes fainter than the LF peak. The mass distribu- 
tion that Paresce & De Marchi (2000) have subsequently 
inferred from this LF shows that the characteristic mass 
of M15 is ~ 0.30 Mq and that the number of stars less 
massive than 0.3 M quickly drops off. This appears to 
be a common feature of all Galactic globular clusters for 
which deep LFs are available, regardless as to their metal- 
licity, position in the Galaxy and dynamics (Piotto et al. 
1997; Paresce & De Marchi 2000). This MF is believed to 
be representative of the IMF of globular clusters, however 
the true IMF can be reliably established only by disentan- 
gling the cluster dynamical evolution from the observed 
MF. This is best achieved when MFs are available at sev- 
eral distances from the cluster centre. For M 15, MFs are 
found in the literature near the core (De Marchi & Paresce 
1996; Sosin & King 1997) and for the above mentioned 
field at 4'.6 NW of the cluster centre. Deep images taken 
by HST/NICMOS Camera 3 (NIC3) during the 1998 par- 
allel campaign have allowed us to derive the MF in three, 
outer fields at 7' NE of the centre. We have used this mass 
distribution together with those published previously to 
constrain the Global Mass Function (GMF) of M 15. If 
the interaction of M 15 with the Galactic tidal field has 
been as weak as recent works indicate (Gnedin & Ostriker 
1997; Dauphole et al. 1996), then the GMF should reflect 
the IMF. 

The NIC3 observations are presented in Section 2 and 
the data reduction is described in Section 3. The LF of 
the NIC3 fields in the H band and its corresponding MF 



are discussed in Section 4 and compared with other liter- 
ature measurements in Section 5. We derive the GMF in 
Section 6 and our conclusions follow in Section 7. 

2. Observations 

M 15 has been observed with the NIC3 camera of 
HST/NICMOS on 1998 July 7th and 18th, during the par- 
allel observations campaign. Three overlapping fields have 
been imaged at about 7' NE from the centre of the cluster, 
at a distance of 7 times the half-light radius (1'; Trager et 
al. 1995). Multiple exposures have been taken of each field 
through both the F110W and F160W filters, centered at 
1.1 /um and 1.6 /im, respectively. The coordinates of the 
fields and the names, filters, and total exposure times of 
the image data-sets are given in Tabic 1. Hereafter we refer 
to the F110W and F160W bands as J and H, respectively. 

3. Data Reduction 

The images have been reduced using the NICMOS stan- 
dard calibration pipeline: they have been first processed 
with CALNICA for bias subtraction, dark-count correc- 
tion and flat-fielding. Images belonging to the same field 
have then been associated by means of the IRAF routine 
MAKEASSOCIATION and combined with CALNICB, to 
remove cosmic rays and to increase the signal-to-noise ra- 
tio. 

Photometry has been performed on each of the three 
combined images with the DAOPHOT package. Stars have 
been identified with DAOFIND, by setting the detection 
threshold at 5cr above the local background. Wc have 
traced the radial profile of each identified object and dis- 
carded those with full width at half maximum (FWHM) 
larger than 2.5 pixels, since the typical FWHM of a well 
defined point source in our frames is 1.5 pixels. Moreover, 
we have compared images of the same association in or- 
der to identify bad pixels not flagged by the calibration 
pipeline. In this way we have selected a sample of 539 stars 
imaged in both J and H bands. Because of the highly 
variable background, we decided to measure stellar count- 
rates in small fixed apertures of 2 pixels in radius (equiva- 
lent to (X'4), and the corresponding background values in 
a fixed annulus with a radius of 5 pixels and a width of 2 
pixels. After background subtraction and before applying 
any aperture correction, we have corrected the count-rates 
for the NIC3 intra-pixel sensitivity, using the equations 
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Table 2. Estimated photometric errors 



Bin 




Error in H 


Error in J 


mag. 


< 22 


0.05 


0.05 


22 < 


mag. < 23 


0.10 


0.09 


mag. 


> 23 


0.12 


0.21 



computed by Storrs et al. (1999; see Table 2) in the case 
of out-of-focus campaign data. 

The aperture correction was determined in three steps: 

i) First, we constructed a mean growth curve for each 
frame from a sample of bright and isolated stars. The 
stellar fluxes were measured in 10 apertures, with radii 
ranging from 1 through to 7 pixels, and the sky was 
taken in a fixed annulus with a radius of 7 pixels and 
a width of 3 pixels. After background subtraction, the 
stellar count-rates obtained for the same aperture were 
averaged into a mean growth curve, from which we 
derived the amount of energy encircled between 2 and 
5 pixels, needed to scale our count rates to an aperture 
of 5 pixels. 

ii) Since NIC3 was out of focus during our observations, 
we used the TinyTim software (Krist & Hook 1999) 
to simulate the instrumental point spread function 
(PSF) with the precise optics settings corresponding 
to a specific filter and observation date. We computed 
two PSFs for each frame, one for our observation date 
(July 1998) and one for 1998 January 15, when NIC3 
was in-focus (in-focus campaigns were carried out in 
January and June 1998). We calculated the encircled 
energy for a 5 pixels aperture for each PSF and used 
the flux ratio of in-focus and out-of-focus to correct 
our measured count rates. 

iii) We finally multiplied the sample count rates 
by the factor of 1.075 so as to correct them 
to the values measured in a nominal infinite 
aperture (NICMOS Photometry Cookbook, cf. 
http: / / www.stsci.edu/hst / nicmos / performance ) . 

The corrected count rates c were then converted to 
magnitudes in the VEGAMAG photometric system by 
means of the relation: 

m = -2.5Zo 9 , (1) 

where U is the conversion factor from flux to count 
rate and Z is the flux for a zero magnitude star in the 
VEGAMAG system, provided for all NICMOS filters 
and VEGAMAG bands by the HST Data Handbook 
( http: / / www.stsci.edu/hst / nicmos/documents/handbooks ) . 

3.1. Photometric uncertainty 

We estimated the internal uncertainty of our photome- 
try by comparing the resulting magnitudes of those stars 
in common to two different fields: Fields 1 and 2 overlap 
nearly over two quadrants sharing 83 stars, whilst Fields 1 



Table 3. Completeness factors 



Field 


bin 


/ (J band) 


/ (H band) 


1 


15-22 


1. 


1. 


2 




1. 


1. 


3 




1. 


1. 


1 


22-23 


0.959 


0.808 


2 




0.949 


0.935 


3 




0.913 


0.864 


1 


23-23.5 


0.807 


0.526 


2 




0.935 


0.821 


3 




0.909 


0.680 


1 


23.5-24 


0.667 


0.274 


2 




0.898 


0.545 


3 




0.870 


0.423 


1 


24-24.5 


0.304 


0.126 


2 




0.640 


0.269 


3 




0.433 


0.118 


1 


24.5-25 


0.132 


0.073 


2 




0.488 


0.122 


3 




0.338 


0.076 



and 3 have in common only 12 stars. As a measure of 
the uncertainty, we used the difference between the mag- 
nitudes measured in each field and the weighted mean of 
the two values, with the weight given by the square root 
of the exposure time. The resulting errors, for Field 2, are 
shown in Table 2 for three magnitude ranges. We associate 
an uncertainty of 0.05 mag to stars brighter than magni- 
tude 22 and a photometric error > 0.1 to fainter objects. 
As regards Field 1, errors are usually twice as large, since 
the exposure time is ~ 3.5 times shorter in both bands. 
For Field 3 we adopted the same uncertainty scale as for 
Field 2. 

3.2. Photometric completeness 

We used the ADDSTAR routine in DAOPHOT to deter- 
mine the completeness of our photometry. We tested each 
frame separately in both the J and H band, by adding 
about 10% of the total number of detected stars in or- 
der not to increase the crowding in the images. We per- 
formed four runs for several magnitude bins. The results 
are shown in Table 3 for all the observed fields. Since NIC3 
is less sensitive in the H band, the completeness in H sig- 
nificantly affects our subsequent data analysis. In the case 
of Field 1, for which only short exposures are available, 
the completeness rapidly decreases to ~ 50 % at H ~ 23. 
Fields 2 and 3, on the other hand, reach a completeness of 
~ 50 % and ~ 40 %, respectively, at H ~ 24. 

4. Analysis of the NICMOS data 

4.1. The colour-magnitude diagram 

We have corrected the observed magnitudes of our sample 
for reddening by assuming E(B-V) = 0.1 from Durrell & 
Harris (1993). This implies Aj = 0.1 and A H = 0.06. 
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Fig. 1. De- reddened CMD of each observed field and all of the fields combined. Stars common to two fields have been 
removed from each field but are included in the CMD of the whole sample, with averaged magnitudes and colours. 



The de-reddened colour-magnitude diagram (CMD) is 
plotted for each of the three fields and for the whole sam- 
ple in Figure 1. For the stars in common to more than one 
field, we have adopted the mean magnitudes and colours 
computed as above. Similar levels of photometric com- 
pleteness are reached for all three fields at magnitudes 
brighter than Ho ~ 23, so that their CMDs can be directly 
compared to one another. In all the three cases the main 
sequence is well defined down to Ho ~ 22 and spreads out 
for —0.4 < (J — H)o < 1.8 at fainter magnitudes due 
to our photometric uncertainty. The few stars brighter 
than H = 18 are probably foreground objects. Indeed, 
Durrell & Harris (1993) determined the turn-off for M 15 
at V = 19.4 which corresponds to H — 18.4 from the 
theoretical track of Baraffe ct al. (1997) at the metallicity 
of M15 ([Fe/H] = -2.15). This implies that the stars at 
Ho — 18 are either cluster objects evolved off the main 
sequence or simply foreground stars. 



The bright tip of the main sequence shows up at Ho — 
18 for Fields 1 and 2 while it falls at H ~ 19 for Field 3. 
This apparently fainter turn-off magnitude is likely due to 
statistical fluctuations in our small sample. 

In order to reduce the contamination from foreground 
and background stars, we applied to the CMD of each 
field a 2.5c clipping selection around the average colour 
of the main sequence. The resulting three decontaminated 
samples were merged onto the CMD of Figure 2 where the 
stars in common are represented with averaged magni- 
tudes and colours. Using the predictions of Ratnatunga 
& Bahcall (1985), we have estimated the field-star con- 
tamination to be about 4 stars for each NICMOS field in 
the direction of M 15. This estimate is valid for a limiting 
magnitude of H = 23 (i.e. 0.2 M in Figure 2), which cor- 
responds to V ~ 27 in the evolutionary tracks of Baraffe et 
al. (1997), and is integrated over the whole (B—V) colour 
range taken into account by Ratnatunga & Bahcall (1985). 
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Therefore, we do not expect the CMD to change signifi- 
cantly after the er-clipping is applied. The photometric er- 
rors are also indicated on the left-hand side of Figure 2: in 
the range 14 < H < 22 the (J — H) colours are known 
with an accuracy of ±0.07 mag, at 22 < H < 23 they 
increase to ±0.13 mag and for Hq > 23 the photometric 
accuracy is as poor as ±0.24 mag. 

Superposed on the observed CMD distribution is the 
theoretical track obtained from the models of Baraffe et 
al. (1997) for a metallicity of [M/H] = -2, which closely 
matches the [Fe/H] = -2.15 value of M 15. This track is 
scaled by the distance modulus of 15.11 mag (Durrell & 
Harris 1993). The stellar masses actually defining the the- 
oretical track are listed on the right-hand side of Figure 2 
for decreasing Hq magnitudes. The main sequence spans 
a mass range between 0.8 M Q and 0.2 M . 

4.2. The luminosity function 

The LF observed for the external fields of M 15 is plot- 
ted in Figure 3 (solid line) corrected for incompleteness. 
The stars in the CMD were grouped into magnitude bins 
(each 0.5 mag wide) between Hq = 18 and Hq — 24, where 
completeness drops to less than 50 %. Since the complete- 
ness of the three fields is nearly the same at Hq < 23, we 
have assumed a mean completeness factor with which we 
have corrected the star counts of bins brighter than 23. 
At Hq > 23, Fields 2 and 3 are the major contributors to 
the observed LF, since they are deeper. In this magnitude 
range we have, therefore, ignored the stars in Field 1 not 
contained in Fields 2 or 3 and have rescaled the total num- 
ber of objects in these two latter fields to match the total 
area of the survey. Therefore we have computed a mean 
completeness factor between these two fields and used the 
resulting value to correct the star counts of bins fainter 
than Hq = 23. Given the shallower photometric depth 
of Field 1 , having ignored its contribution to the LF for 
Hq > 23 is not likely to affect the statistical significance 
of our results. 

The LF of Figure 3 extends over the range 2 < Mh < 
9, peaking at Mh = 7.1 (corresponding to Mj = 8.5; 
Baraffe et al. 1997) and is characterised by a quite sharp 
drop to fainter magnitudes due to a lack of progressively 
less massive stars. Moreover, this LF confirms the gen- 
eral behaviour found by Paresce & De Marchi (2000) for 
a dozen Galactic globular clusters: their average LF, com- 
puted from stars below 1 M Q near the half-light radius, 
rises to a maximum value at Mj ~ 8.5 — 9 and then drops 
for fainter magnitudes (this applies regardless of the clus- 
ter position and orbit in the Galaxy and of its internal 
dynamical state). 

4.3. The mass function 

The MF of the observed fields was derived from the LF of 
Figure 3. Instead of deriving the MF by inverting the LF, 
we followed the approach of Paresce & De Marchi (2000), 




0.7 
0.6 
0.5 

0.4 
0.3 



(J H)„ 



Fig. 2. De-reddened CMD of the whole sample on which 
the theoretical track of Baraffe et al. (1997), computed for 
the metallicity of M 15 ([Fe/H] = —2.15), has been super- 
posed. Photometric errors are represented, as a function of 
Ho, at the left-hand side of the diagram, whilst the mass 
range spanned by main sequence stars is marked on the 
right-hand side of the plot, in unit of solar masses. 

so as to treat separately observational and theoretical un- 
certainties. We assumed a model MF of log-normal type, 
i.e. one of the type: 



«( 



dN 



\dlog(m) 



= A 



log(m/m c 



1 2 



V2(T 



(2) 



with characteristic mass m c and standard deviation a. A is 
a normalization constant. We then folded it through the 
derivative of the mass-luminosity relationship of Baraffe 
et al. (1997) to obtain a model LF, which we compared 
to the data until a suitable value of the parameters was 
found that gives a good fit to the observations. The best 
fitting LF is superposed to the observations in Figure 3 as 
a dashed distribution and is obtained with m c = 0.3 M© 
and a = 0.28. 

5. A comparison with existing surveys of M 15 

As mentioned in the Introduction, M 15 has been exten- 
sively studied for it is at a relatively large distance from 
the Galactic plane and, consequently, not severely affected 
by dynamical interactions with the Galaxy. 

Figure 4 spatially visualises the stellar fields observed 
in the core of M 15 by De Marchi & Paresce (1996; 
DMP96) and Sosin & King (1997; SK97-1, SK97-2) in 
order to derive the central MF. The DMP96 and SK97-1 
fields are at the centre, whilst the SK97-2 field is 20" away 
from it. All fields are well within the half-light radius of 



6 



A. Pasquali et al.: The global mass function in M 15 



100 



a 
'2 

ab 

E 




Fig. 3. Luminosity function (solid curve) of the 
whole sample, corrected for photometric completeness. 
Poissonan errors have been associated to the observed 
star counts. The dashed distribution is the best fitting 
LF, computed under the assumption of a log-normal mass 
distribution with a characteristic mass of 0.3 M©. 
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Fig. 5. Luminosity functions as derived by De Marchi 
& Paresce (1996; solid curve) and Sosin & King (1997; 
dashed distribution) for their field positioned at the cen- 
tre of M 15. Poisson errors have been associated with the 
plotted star counts. 



Table 4. M 15 surveys 




Fig. 4. Spatial distribution of the stellar fields observed 
within the half-light radius of M 15 (l';Trager et al. 1995). 
DMP96 is the field studied by De Marchi & Paresce (1996) 
near the core, while SK97-1 and SK97-2 are the areas sur- 
veyed by Sosin & King (1997). SK97-2 is 20" away from 
the centre. 



M15 (r h ~ 1'; Tragcr et al. 1995). Both studies detect 
a substantial amount of mass segregation, although the 
LFs differ for My > 4.8. The two luminosity distributions 
are plotted in Figure 5: the LF of De Marchi & Paresce 
(1996) has been here translated from the original FOC 





Distance from 


Filters 




the center 




Durrcll & Harris (1993) 


7' NW 


B,V 


De Marchi & Paresce (1995) 


4'.6 NW 


v,i 


this work 


7' NE 


J,H 



F346M band magnitudes to standard Johnson V values 
using Baraffc et al.'s tracks (private communication). 

Both LFs have been normalised to the peak star count. 
An interesting point is that at My > 4.8 the LF of De 
Marchi & Paresce (1996) dramatically drops, whilst the 
LF derived by Sosin & King (1997) is flat. We believe that 
this discrepancy is due to photometric incompleteness and 
to the colour extrapolation from the UV to the V band. 

Whilst mass segregation affects the innermost region of 
globular clusters, dynamical interactions with the Galaxy 
act predominantly on their outermost regions inducing 
stellar evaporation and stripping. Thus, clusters as rich 
as M 15 are likely to be dynamically unperturbed at their 
half-light radius (Richer et al. 1991) and stars at this dis- 
tance can be used to constrain the IMF. Table 4 lists the 
observations available in the literature which cover the 
outskirts of M 15. Their corresponding fields are plotted 
in Figure 6 together with the cluster centre and half-light 
radius. 

We extracted the LFs from these papers and ho- 
mogenised them in units of I magnitudes by using Baraffc 
et al.'s (1997) tracks. In Figure 7 we over-plot the LF de- 
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Fig. 6. Spatial distribution of the fields observed at dis- 
tance larger than the cluster half-light radius. DH93 labels 
the field imaged by Durrcll & Harris (1993) at 7' NW, 
whilst DMP95 represents the field observed by De Marchi 
& Paresce (1995) at 4.'6 NW from the centre, fl, f2 and £3 
are the NIC3 fields analysed in this paper. 

rived from the NIC3 data on Durrell & Harris' (1993; 
left-hand panel) and De Marchi & Paresce's (1995; right- 
hand panel) distributions. The latter were scaled to the 
peak star counts of the NIC3 distribution. Unfortunately, 
Durrell & Harris' (1993) observations are not deep enough 
to reach the LF turn-over at Mi ~ 8.5. Nevertheless, their 
LF agrees well with the ascending part of the luminos- 
ity distribution derived from the NIC3 exposures. Very 
good overlap is instead achieved between the LFs of this 
paper and that of De Marchi & Paresce (1995) over the 
range 4 < Mj < 10, so that both resulting MFs peak near 
0.3M Q . 

6. Dynamical structure 

Having derived the spatially resolved LF for M 15, we can 
now study its radial changes and address the issue as to 
whether they are consistent with mass segregation ensuing 
from two body relaxation. To study the dynamical prop- 
erties of the cluster, we employed the multi-mass Michie- 
King models originally developed by Meylan (1987, 1988) 
and later suitably modified by Pulone et al. (1999) and De 
Marchi et al. (2000) for the general case of clusters with 
a set of radially varying LFs. Each model run is charac- 
terised by a MF in the form of a power-law dN / d log m oc 
m~ x , with a variable index x, and by four structural pa- 
rameters describing, respectively, the scale radius (r c ), the 
scale velocity (v s ), the central value of the dimensionless 
gravitational potential W , and the anisotropy radius (r a ). 
(After having suggested, in Section 3, a rather general de- 
scription of the functional form of the MF, namely a log- 
normal distribution, it might seem inappropriate to adopt 
a variable power-law as the basis for the MF in our dy- 



namical model. In fact, as we show below, the resulting 
MF is indistinguishable, for any practical purposes, from 
a log-normal distribution.) 

From the parameter space defined in this way, we se- 
lected those models that simultaneously fit both the ob- 
served surface brightness profile (SBP) and velocity dis- 
persion profile (VDP) of the cluster as measured, respec- 
tively, by Guhathakurta et al. (1996; for r < 100") and 
Trager et al. (1995; for r > 100") and by Gebhardt et 
al. (2000). However, even requiring good fits to both the 
SBP and VDP can, by itself, only constrain r c , v s , W , 
and r a , whilst still allowing the MF to take on a variety 
of shapes. To break this degeneracy, we further imposed 
the condition that the model MF agree with the observed 
LF at all radial distances offered by the data. 

Since Michie-King modeling only provides a "snap- 
shot" of the current dynamical state of the cluster, it is 
useful to refer to the GMF, i.e. the mass distribution of all 
cluster stars at present, or, in other words, the MF that the 
cluster would have simply as a result of stellar evolution 
(i.e. ignoring any local modifications induced by internal 
dynamics and/or the interaction with the Galactic tidal 
field). Clearly, in this case the IMF and GMF of main 
sequence (unevolved) stars is the same. For practical pur- 
poses, the GMF has been divided into sixteen different 
mass classes, covering main sequence stars, white dwarfs 
and heavy remnants, precisely as described in Pulone et 
al. (1999). 

Our parametric modelling approach assumes energy 
equipartition amongst stars of different masses. Thus, we 
ran a large number of trials to see whether we could find 
a set of parameters for the GMF (i.e. a suitable GMF 
"shape") such that the local MFs produced by mass seg- 
regation would locally fit the observations. We note here 
that, rather than converting the observed LFs into MFs for 
comparison with the predictions of the model, we prefer 
to keep observational errors and theoretical uncertainties 
separate. Therefore, we convert the model MFs to LFs us- 
ing for all the same M-L relation, namely that of Baraffe 
et al. (1997), precisely as we did in Figure 3. Not surpris- 
ingly, our exercise confirms what we had already shown 
in that Figure and described above: as long as a single 
value of the index x is used for the GMF over the mass 
range 0.2 — 0.8 M Q , none of the predicted local LFs can be 
fitted to our data. In fact, a change of slope is needed at 
m ~ 0.3 Mq so that both the rising and dropping portions 
of the local LF can be simultaneously reproduced. If we 
then allow the MF to take on more than one slope, the 
GMF that best fits the observations is one with x = 0.9 
for stars in the range 0.3 — 0.8 M Q and x = —2.2 at smaller 
masses. The shape of this GMF is, thus, very similar to 
the log-normal distribution shown in Figure 3. 

The set of LFs predicted by the set of Michie-King pa- 
rameters that best fit all available observations is shown in 
Figure 8, where the squares correspond to the LFs avail- 
able in the literature for this cluster at various distances 
from its centre. The fit to the SBP and VDP obtained with 
the same set of parameter values is shown in Figures 9 
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Fig. 7. The LFs computed by Durrcll & Harris (1993, DH93) and De Marchi & Paresce (1995, DMP95) are compared 
with the distribution obtained from the NIC3 data. All of the luminosity distributions have been translated into Mj 
values and scaled to the peak star-counts of the NIC3 luminosity function. 



and 10 and is surprisingly good. The values of the best 
fitting structural parameters are shown in Table 5, where 
they can be compared with those in the literature. The 
agreement is excellent, apart from a small difference in 
the value of the tidal radius which is, admittedly, not seri- 
ously constrained by our data. We note here that we can 
directly compare the observed SBP with our model since 
the solid line in Figure 9 corresponds to stars of ~ 0.8 M Q , 
namely those contributing most of the cluster's light. As 
one should expect, stars in different mass classes have dif- 
ferent projected radial distributions. 

Although stars more massive than ~ 0.8 M Q have 
evolved and are no longer visible, the shape of the IMF 
in this mass range has strong implications on the fraction 
of heavy remnants in the cluster and, as such, on the cen- 
tral velocity dispersion. We find that the best fit to the 
data and to the cluster's structural parameters, as given 
above, requires a value of x = 0.8 for stars in the range 
100 — 0.8 M . It should be noted that the global cluster 
MF is thus slightly shallower than Salpeter's IMF, which 
would have x = 1.35. The total implied cluster mass is 
<~ 5.4 x 10 5 M Q and the mass-to-light ratio is on average 
M/L ~ 2.1, with M/L ~ 3.7 in the core. The total cluster 
luminosity Ly has been estimated by integrating the best- 
fitting SBP (solid line in Figure 9), properly normalised to 
match the observed central surface brightness. The best 
fitting models suggest that a large fraction of mass (of 
order <~ 44%) is trapped in heavy remnants, namely stel- 
lar black holes, neutron stars and white dwarfs. However, 
the presence of an intermediate mass black hole is not 
required. 

Interestingly, the rather shallow GMF that we obtain 
for stars above 0.8 M , which, in turn, results in a large 
fraction of heavy remnants, is also dictated by the cen- 
tral enhancement seen in the SBP. The excellent fit that 
our model offers to the radial surface density of TO-mass 
stars (Figure 9) would not be otherwise possible. A steeper 



Table 5. Parameters of the Michie - King models used 
for M 15 



Parameter 


Fitted 


Literature 


Ref. 




value 


value 




core radius r c 


2"3 


2" 


a 


tidal radius r t 


17' 


2i:5 


b 


half-light radius r h 


lf3 


1' 


c 


central vel. disp. ov 


12.1 km s _1 


11.7 km s _1 


d 



a: Guhathakurta et al. (1996) 

b: Harris (1996) 

c: Djorgovski (1993) 

d: Gebhardt (2000) 



MF index would result in a shallower central profile which 
would fail to reproduce the central density enhancement. 
It should also be noted that, since a canonical King- type 
profile does not reproduce the SBP of M15 (see, e.g., 
Trager et al. 1995), the value of r c that we obtain does 
not correspond with the canonical definition of core radius, 
namely that at which the surface density falls to one-half 
its central value. The latter definition only applies to the 
profile of the original King model with a single mass com- 
ponent (see Kent & Gunn 1982). In our case, r c simply 
represents a scale radius. 

7. Other Michie-King models for M15 

Contrary to what we have concluded here, in their analysis 
of the radial variation of the LF of M15, Sosin & King 
(1997) concluded that a multi-mass Michie-King model is 
unable to reproduced the observations. As their Figure 12 
shows, the variation predicted by their multi-mass model 
is larger than that observed when comparing the centre of 
the cluster and the region at r ~ 5'. We have identified 
three reasons that might have led Sosin & King (1997) to 
this conclusion and we discuss them here briefly. 
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Fig. 8. The best fits of the Michie-King models to the LF observed at 7' (this work), 4.' 6 (De Marchi & Paresce 1995), 
and 20" (Sosin & King 1997) from the cluster centre and in the core (De Marchi & Paresce 1996; Sosin & King 1997). 



The first is mostly related to the approach used and, 
as such, should only affect the uncertainty of the results. 
Sosin & King (1997) decided to transform the observed 
LFs into MFs, not necessarily using the same M-L rela- 
tion for all data, and to compare the predictions of their 
multi-mass models to these MFs. The advantage of our 
approach, in which the LF predicted by the model is com- 
pared with the observed LF, is that we ensure that obser- 
vational errors and theoretical uncertainty (in the model 
and M-L relation) are kept separate and that only one 
M-L relation is used throughout the process. 

Secondly, it appears that their dynamical model is un- 
able to reproduce at all the observed velocity dispersion 
profile, as the authors themselves point out. Conversely, 
ours is in excellent agreement with the observations. As 
a result, Sosin & King (1997) predict a fraction of heavy 
remnants (< 1%) well below the current estimates for this 
cluster. 

Most importantly, however, we believe that their in- 
ability to reproduce the observed radial variation of the 
LF stems from the functional form of the MF that they 
adopt. As we describe in Section 4, we have made a general 
assumption about the shape of the GMF, in the form of 



a log-normal distribution, based on what was learnt from 
the observation of a large number of halo GCs (Paresce & 
De Marchi 2000). We then let our procedure find the pa- 
rameter values that simultaneously fit all available data. 
Since the number of independent measurements is larger 
than that of the unknowns, the procedure is bound to con- 
verge. On the other hand, Sosin & King (1997) adopt the 
MF determined by Piotto ct al. (1996) at r ~ 5' as the 
basis for their dynamical model, but the MF predicted by 
it for the cluster core fails to match the data. The origin of 
the mismatch seems to lie predominantly in the assumed 
shape of the MF, which is flat in the range m > 0.7 M 
and then sharply rises at lower masses. This results in a 
MF in the central cluster regions that sharply drops in 
the range 0.8 M > m > 0.7 M Q . It appears that having 
adopted a model MF with a more gentle rise from 0.8 M Q 
all the way through to 0.5 M Q (which still fits the data at 
5' equally well) would have produced a MF in consider- 
ably better agreement with the observations in the central 
cluster region. 
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Fig. 9. Model fit to the surface brightness profile. The solid line corresponds to the profile of stars of ~ 0.8 M Q , 
responsible for most of the light of the cluster. 



8. Conclusions 

We have analysed HST/NICMOS parallel data available 
for three fields in the Galactic globular cluster M 15 lo- 
cated at 7' NE of the cluster centre. Their total LF ex- 
tends over the range 2 < Mh < 9 and peaks at Mh = 7.1 
(or Mi ~ 8.5). It also shows a sharp drop towards fainter 
magnitudes, which is a typical signature of the lack of pro- 
gressively less massive stars. We have fitted this LF with 
a log-normal mass distribution and obtained a character- 
istic mass of 0.3 M Q , with a standard deviation of 0.28. 
These values agree well with the characteristic mass and 
standard deviation derived for a dozen Galactic globular 
clusters by Paresce & De Marchi (2000), thus support- 
ing the ubiquity of the log-normal mass distribution for 
globular clusters. 

The H-band LF obtained at 7' from the centre of M 15 
has been compared with the LFs derived by Durrell & 
Harris (1993) and De Marchi & Paresce (1995) at 7' NW 
and 4.'6 NW of the cluster centre, respectively. The com- 
parison has required the translation of all original ob- 
served magnitudes into the I band. The overlap among 
these three LFs is excellent, indicating that the distribu- 



tion of the stars at distances larger than the half-light 
radius (~ 1'; Trager et al. 1995) may not be significantly 
perturbed in M 15, as expected from its Galactocentric 
distance (Rq — llKpc; Gnedin & Ostrikcr 1997) and 
the small ellipticity of its orbit (e = 0.32; Dinescu et al. 
1999). Indeed, calculations by Gnedin & Ostrikcr (1997) 
show that the time to destruction of M 15, due to the com- 
bined effects of internal dynamical evolution and interac- 
tion with the Galaxy, is as large as 50Gyr. Therefore, if 
correct, this time-scale would suggest that the outermost 
fields observed in M 15 have not been significantly per- 
turbed by tidal stripping and evaporation and that their 
content should likely represent the initial stellar mass dis- 
tribution. In other words, the high degree of similarity 
between the LFs (and hence the MFs) of the fields at 4.' 6 
and 7' from the centre would imply that these are very 
close to be the cluster IMF and there are no significant 
radial variations in the cluster IMF. 

We have used the LF measured for M 15 at 7' (this 
work), 4.'6 (De Marchi & Paresce 1995), and 20" (Sosin 
& King 1997) from the cluster centre and in the core (De 
Marchi & Paresce 1996; Sosin & King 1997) to study the 
effects of mass segregation. We have fitted Michie-King 
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Fig. 10. Model fit to the velocity dispersion profile. 



models to the observed surface brightness and velocity 
dispersion profiles in order to estimate the cluster struc- 
tural parameters and to the observed LF to constrain the 
shape of the cluster GMF. The latter turns out to be char- 
acterised by two slopes, x = 0.9 for stars in the range 
0.3 — 0.8 M Q and x = —2.2 at smaller masses, and is thus 
very close to the log-normal distribution obtained directly 
from our NIC3 data near the cluster's half-light radius. 

The values of the cluster structural parameters that 
best fit the observations imply a cluster total mass of 
— 5.4 x 10 5 M Q and a mass-to-light ratio of M/L ~ 2.1 on 
average, with M/L ~ 3.7 in the cluster core. In addition, 
the best-fitting Michie-King model parameters suggest a 
slope of x = 0.8 for the IMF in the range 100 - 0.8 M Q , 
which supports the presence of a large fraction of heavy 
remnants (~ 44%) in the core. If such a high fraction 
of heavy remnants is present, as originally suggested by 
Illingworth & King (1977), it would rule out the need of 
an intermediate-mass black hole to explain the mass seg- 
regation and velocity dispersions observed in the core of 
M15. 
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